Vorticity generation at second order in cosmological perturbation theory 
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We show that at second order in cosmological perturbation theory vorticity generation is sourced 
by entropy gradients. This is an extension of Crocco's theorem to a cosmological setting. 
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I. INTRODUCTION 

It is well known in fluid dynamics that vorticity gen- 
eration is sourced by entropy gradients. This was first 
pointed out by Crocco in 1937 Despite vorticity be- 
ing ubiquitous in nature, studies of vorticity in the early 
universe and in cosmology have so far been rare. 

At linear order in perturbation theory scalar and vec- 
tor perturbations, classified according to their transfor- 
mation behaviour on spatial 3-hypersurfaccs, decouple 
from each other. Scalar perturbations are much easier 
to treat mathematically and play the dominant role in 
structure formation on super-horizon scales. At linear 
order vorticity, intrinsically of vectorial nature (in fluid 
dynamics simply the curl of the fluid velocity) cannot 
be constructed from scalar quantities, and the vorticity 
tensor constructed from vector perturbations is, in gen- 
eral, sourced only by anisotropic stress and decays in its 
absence @,i,i,i,[6|. ^ 

Things are different though at second order in the 
perturbations. Only recently, with second order cos- 
mological perturbation theory becoming mature (see 
e.g. Refs. @, i, 0, [HI, [H, [ll, [3), has the issue of 
vorticity generation beyond the standard linear order 
begun to be addressed [l^, [H, [l^- However, these 
previous studies have been restricted to barotropic fluids 
with an equation of state P = P{p). We show that 
allowing for a non-adiabatic pressure perturbation gives 
qualitatively different results. 

We focus here on the case of vorticity generation 
at second order in the perturbations, sourced only by 
scalar and vector perturbations. We consider a perfect 
fluid, namely a fluid whose energy-momentum tensor is 
diagonal (the inclusion of anisotropic stress will source 
vorticity already at first order in the perturbations). 
Such a fluid has an equation of state P = P{p,S), 
where P is the pressure, p is the energy density and 
S is the entropy of the fluid. We can expand the 
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pressure perturbation for a fluid with this equation of 
state as (5P = If (JS* 



SP = <5P„ad + ciSp, 
speed, defined as Cg e 



f ^Sp, which can be written as 

where Cg is the adiabatic sound 
E P'/p', and we have defined the 
non-adiabatic pressure perturbation i5Pnad = ff Ip'^'^'j 
in addition to the pressure and energy density pertur- 
bations. For a detailed discussion of the non-adiabatic 
pressure perturbation, see Ref. [3 • 

Our main result, derived in detail below and given in 
Eq. (|3.6|1 . can be written concisely as 



^Lj K Spi ySPni^dl, 



(1.1) 



that is, the gradients in the non-adiabatic pressure per- 
turbation coupled to gradients in the density act as a 
source for vorticity at second order. 

We use cosmological perturbation theory throughout, 
which will make the numerical implementation of the re- 
sults straightforward and, compared to other approaches, 
has the added benefit of actually being physically trans- 
parent. The Paper is organised as follows. In the 
next section we define and introduce the key variables. 
In Section IIIII we define the vorticity tensor and give 
its evolution. We discuss our results and conclude in 
Section IIVI All the governing equations necessary to 
derive the results in Section|TTT]are given in the appendix. 

In this Paper, we use conformal time, rj, throughout, 
denoting derivatives with respect to conformal time with 
a prime. The scale factor is a, and the Hubble parameter 
is 7i, where Ti. = a' /a. Greek indices, pi,v,\ run from 
0, . . . , 3, while lower case Latin indices, i, j, fc, take the 
value 1,2, or 3. Covariant derivatives are denoted by a 
semi-colon, partial derivatives by a comma. The order 
of the perturbations is denoted with a subscript immedi- 
ately after a perturbed quantity. We work in the uniform 
curvature gauge throughout. 



II. DEFINITIONS 



^ Note, however, that there are certain situations in which vorticity 
is not sourced by anisotropic stress at linear order. See Ref. [3 
for an example of first order vorticity sourced by heat flux. 



The Friedmann-Robertson- Walker (FRW) metric ten- 
sor, up to and including second order perturbations, has 
in the uniform curvature gauge the components (see. 
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e.g. Ref. [13) 
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-a^ (1 + 201 + 02) , 



a' I Bii + 2^2! 



(2.1) 



where we have assumed a flat {K = 0) background with- 
out loss of generality. We neglect tensor perturbations, 
which will add another source term to the momentum 
conservation equation, but are beyond the scope of this 
work [l9[. Here a = a{r]) is the scale factor, 01 and 02 are 
the lapse functions at first and second order, respectively, 
and Bii and B2i, represent the shear in this gauge. All 
perturbed quantities are function of x^. Note Bu and 
B2i can be further split into scalar and divergence-free 
vector parts p^ . though this step is unnecessary here. 

The fluid four velocity, u'^, obeying the normalisation 
condition m^m^ = — 1, has components 
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Vu + 2^2^ - 01 



, (2.2) 
(2.3) 



up to second order, where is the spatial three velocity, 
and is the covariant spatial velocity, defined as = 



III. VORTICITY 

The vorticity tensor is defined as the projected anti- 
symmetrised covariant derivative of the fluid four veloc- 
ity, that is 01 



(3.1) 



where U[a;/3] = 5 (uq;/3 — u/3-a), and the projection tensor 
V^u into the instantaneous fluid rest space is given by 



(3.2) 



The vorticity can then be decomposed, up to second 
order in perturbation theory, as ojij = ujnj + ^uj2ij , where 
the first order part is simply 



and the second order part is 



(3.3) 



uj2ij = a^2[j,i]+2a K[j^ij] + 01, {Vi + B^)^^ - 0i5i[jj] 

(3.4) 

Using the evolution equations given in Appendix [Aj 
the first order part evolves as 



SHc^ujuj = , 



(3.5) 



which gives the well known result that at first order, 
in the absence of an anisotropic pressure source term, 
\ujiijujj^''\ cx during radiation domination, where 

Cg = 1/3 (3]. Hence the vorticity remains zero in this 
case, if it is initially zero. 

However, at second order we get a non-zero source 
term for the vorticity evolution equation, assuming zero 
anisotropic pressure. ^ Even assuming zero first order vor- 
ticity, ujiij = 0, the second order vorticity evolves accord- 
ing to 
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3'Hc'^uj2ij 

2a r 
' po + Po{ 



(3.6) 



nad 1 



'5Pl,[j^-Pnadl,j] 

Pq + Pq 



Thus, we sec that the non-adiabatic pressure perturba- 
tion gradients coupled to density perturbation gradients 
act as a source for second order vorticity. In the case 
of a vanishing entropy perturbation, as is the case for 
barotropic fluids, we recover the result of Ref. [l3l. For 
completeness, we give the full second order vorticity evo- 
lution equation without assuming touj = in the Ap- 
pendix as Eq. 



IV. DISCUSSION AND CONCLUSIONS 

In this work we have shown that at second order in 
the perturbations vorticity is generated from first order 
scalar and vector perturbations for a perfect fluid. This 
is an extension of Crocco's Theorem to an expanding, 
dynamical background, namely a FRW universe. 

Whereas in previous works barotropic flow was as- 
sumed, allowing for entropy gives a qualitatively novel re- 
sult. This implies that the description of the cosmic fluid 
as a potential flow, which works exceptionally well at first 
order in the perturbations, will break down at second 
order for non-barotropic flows. Similarly in barotropic 
flow Kelvin's theorem guarantees conservation of vortic- 
ity. This is no longer the case if the flow non-barotropic. 

Entropy perturbations will arise in many settings, 
here we only considered the case of a single fluid 
with non-zero intrinsic entropy. Even in simple single 
field inflation models, after the end of the slow-roll 
regime, there is a non-zero entropy perturbation (or 
non-adiabatic pressure; see, e.g., Ref. [III). Further- 
more, in the case of multiple fluids we get in addition 
to the intrinsic entropies of the individual fluids also 
an entropy component stemming from the mixing of 
the fluids (the relative entropy perturbation in the 
terminology of Ref. [T3|). even if there is no energy and 
momentum transfer between the fluids. This relative 



^ The calculation of the vorticity evolution equation uses first 
order evolution and field equations which we give in Appendix 
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entropy perturbation can be another source of vorticity 
in the multi-fluid case. However in the cosmic plasma it 
is likely that the relative entropy perturbation between 
fluids is subdominant compared to those mechanisms 
discussed below. 

Our result, that the vorticity is non-zero at second or- 
der in the presence of non-adiabatic or entropy pertur- 
bations, has immediate implications for the generation 
of magnetic fields in the early universe, since Biermann 
showed that the generation of magnetic fields is related 
to vorticity [13] (see also Harrison, Ref. [llj). 

Previous works either used momentum exchange 
between multiple fluids to generate vorticity, as in 
Refs. [H m, [23, H Hi, S lii, or used intermediate 
steps to first generate vorticity for example by using 
shock fronts as in Ref. [13]. The vorticity generated in 
the latter case was then used to source turbulence. As 
we have shown above, we do not require such extra steps. 
We will study the evolution of magnetic fields in the case 
of non-zero entropy perturbations, and including tensor 
perturbations, numerically in a forthcoming paper [l9l |. 

At first glance, one might expect that because the vor- 
ticity is generated as a second order effect by products 
of first order density and entropy perturbations, the vor- 
ticity will be of order lO^^*'. However, the source term 
is in fact constructed from gradients of the density and 
entropy perturbations and. because of this, the effect can 
readily be larger than 0(10""'^°). In Fourier space, this 
source term will be 'boosted' by the wavenumbers on 
small scales. Since our mechanism will manifest itself on 
sub-horizon scales and hence the wavenumbers are large, 
this has the potential to increase the effect. For example, 
the authors of Ref. [11] recently showed, in the context 
of second order tensor perturbations, that on sub-horizon 
scales this effect can 'boost' second order quantities to be 
of comparable amplitude as first order Quantities. In fact, 
the WMAP five year results results [20] give only a non- 
zero upper bound on entropy perturbations on horizon 
scales. Since our effect is on sub-horizon scales, following 
the argument of WMAP5 and assuming a blue spectrum 
for the entropy perturbation, it is likely that this effect 
is non-negligible. 

One possible observational signature of vorticity in 
the early universe is that of B-mode polarisation of the 
cosmic microwave background (CMB) radiation. At 



linear order in perturbation theory, only tensor pertur- 
bations (or gravitational waves) will produce B-mode 
polarisation since scalar perturbations only produce 
E- modes (soj . and vector perturbations will become 
subdominant during inflation, and will decay with the 
expansion of the universe. However, at second order the 
former is no longer true and, as we have shown above, 
density perturbations can generate vorticity. This could 
then produce B-mode polarisation large enough to be ob- 
servable by future CMB experiments such as Planck [3l[ . 
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APPENDIX A 

Here we give the evolution and constraint equations 
necessary to derive the results of Section IIIIl For more 
details see Ref. [l^ ]. 



1. Energy- momentum conservation 

In the background the energy conservation equation is 

p'^^-m{p(, + p^) . (Ai) 

At first order we have for the energy conservation equa- 
tion 

5p[+m{5pi+5Pi) + (po + AXfe = 0, (A2) 
and for the momentum conservation equation 



vu + n{l-?,cl)v^, 



Po + Po 



= 0. 



(A3) 



At second order we only need the momentum conserva- 
tion equation 



{pO+Po)V2^ 



+ m (po + Po) + SP2 + (po + Po) h 

-201 



+8H {Spi + SPi) Vu - 2 (po + Po) [^'iBu + {4,1) 

+2 [pa + Po) vl + 2{po + Po) v\Vu,k - 2 (po + Po) vikBl 

I 



+ 2(Spi+SPi)4,,, + 2 {Sp, + 6Pi)Vu (A4) 

i L J 

(po + Po) VI, + (pf, + + 47^ (po + Po) ) Vu 
0. 



2. Field equations 

At zeroth order the evolution of the scale factor is gov- 
erned by the Fricdmann equation, given by 

^""^ ^ (A5) 



We only need the field equations up to first order in 
the perturbations, a considerable simplification, since the 



-a^Po ■ 
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second order lapse function in Eq. (|A4[) above cancels 
when calculating the vorticity evolution. The Einstein 
constraint equations at first order arc 



3. Vorticity 



6W 



and 



(A6) 



V^Bu - fc., - 47^01,, = IGnGa^po + Po)Vu ■ (A7) 



Using the definition of the vorticity tensor, Eq. (|3.4p 
above, and the equations given in the previous subsec- 
tion, we get the evolution equation for the vorticity ten- 
sor at second order 



2ij 



V Po + Po 

+2 {V^^- - ) iuuj.k - 2 {Xl^ - V^^) uuk + 2 {Xl, - l/i^ u^^k 

= r~Fr I 3^ (Vli(5Pnadlj - VijSPnMl,i) H {6Pl jSPnadl.t - Sp^ ,(5Pnadlj) 1 , 

Po + -t^o I Po + ^0 J 



(A8) 



where Xu is given entirely in terms of matter perturba- 
tions as 



Xu = V- 



AttG 



{3H{po + Po)Vu - -5/51^0 



(A9) 
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